home *** CD-ROM | disk | FTP | other *** search
/ Almathera Ten Pack 2: CDPD 1 / Almathera Ten on Ten - Disc 2: CDPD 1.iso / pd / 101-125 / 105 / bison / warshall.c < prev   
C/C++ Source or Header  |  1995-03-13  |  3KB  |  127 lines
  1. /* Generate transitive closure of a matrix,
  2.    copyright (C) 1984 Bob Corbett and Richard Stallman
  3.  
  4.    Permission is granted to anyone to make or distribute verbatim copies of this program
  5.    provided that the copyright notice and this permission notice are preserved;
  6.    and provided that the recipient is not asked to waive or limit his right to
  7.    redistribute copies as permitted by this permission notice;
  8.    and provided that anyone possessing an executable copy
  9.    is granted access to copy the source code, in machine-readable form,
  10.    in some reasonable manner.
  11.  
  12.    Permission is granted to distribute derived works or enhanced versions of
  13.    this program under the above conditions with the additional condition
  14.    that the entire derivative or enhanced work
  15.    must be covered by a permission notice identical to this one.
  16.  
  17.    Anything distributed as part of a package containing portions derived
  18.    from this program, which cannot in current practice perform its function usefully
  19.    in the absense of what was derived directly from this program,
  20.    is to be considered as forming, together with the latter,
  21.    a single work derived from this program,
  22.    which must be entirely covered by a permission notice identical to this one
  23.    in order for distribution of the package to be permitted.
  24.  
  25.  In other words, you are welcome to use, share and improve this program.
  26.  You are forbidden to forbid anyone else to use, share and improve
  27.  what you give them.   Help stamp out software-hoarding!  */
  28.  
  29. #include <stdio.h>
  30. #include "machine.h"
  31.  
  32.  
  33. /* given n by n matrix of bits R, modify its contents
  34.    to be the transive closure of what was given.  */
  35.  
  36.   void
  37. TC(R, n)
  38.   unsigned *R;
  39.   int n;
  40. {
  41.   register int rowsize;
  42.   register unsigned mask;
  43.   register unsigned *rowj;
  44.   register unsigned *rp;
  45.   register unsigned *rend;
  46.   register unsigned *ccol;
  47.  
  48.   unsigned *relend;
  49.   unsigned *cword;
  50.   unsigned *rowi;
  51.  
  52.   rowsize = WORDSIZE(n) * sizeof(unsigned);
  53.   relend = (unsigned *) ((unsigned) R + n*rowsize);
  54.  
  55.   cword = R;
  56.   mask = 1;
  57.   rowi = R;
  58.   while (rowi < relend)
  59.     {
  60.       ccol = cword;
  61.       rowj = R;
  62.  
  63.       while (rowj < relend)
  64.     {
  65.       if (*ccol & mask)
  66.         {
  67.           rp = rowi;
  68.           rend = (unsigned *) ((unsigned) rowj + rowsize);
  69.  
  70.           while (rowj < rend)
  71.         *rowj++ |= *rp++;
  72.         }
  73.       else
  74.         {
  75.           rowj = (unsigned *) ((unsigned) rowj + rowsize);
  76.         }
  77.  
  78.       ccol = (unsigned *) ((unsigned) ccol + rowsize);
  79.     }
  80.  
  81.       mask <<= 1;
  82.       if (mask == 0)
  83.     {
  84.       mask = 1;
  85.       cword++;
  86.     }
  87.  
  88.       rowi = (unsigned *) ((unsigned) rowi + rowsize);
  89.     }
  90. }
  91.  
  92.  
  93. /* Reflexive Transitive Closure.  Same as TC
  94.    and then set all the bits on the diagonal of R.  */
  95.  
  96.   void
  97. RTC(R, n)
  98.   unsigned *R;
  99.   int n;
  100. {
  101.   register int rowsize;
  102.   register unsigned mask;
  103.   register unsigned *rp;
  104.   register unsigned *relend;
  105.  
  106.   TC(R, n);
  107.  
  108.   rowsize = WORDSIZE(n) * sizeof(unsigned);
  109.   relend = (unsigned *) ((unsigned) R + n*rowsize);
  110.  
  111.   mask = 1;
  112.   rp = R;
  113.   while (rp < relend)
  114.     {
  115.       *rp |= mask;
  116.  
  117.       mask <<= 1;
  118.       if (mask == 0)
  119.     {
  120.       mask = 1;
  121.       rp++;
  122.     }
  123.  
  124.       rp = (unsigned *) ((unsigned) rp + rowsize);
  125.     }
  126. }
  127.